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The exactly solvable Kitaev honeycomb lattice model is realized as the low energy effect Hamil- 
tonian of a spin-1/2 model with spin rotation and time-reversal symmetry. The mapping to low 
energy effective Hamiltonian is exact, without truncation errors in traditional perturbation series 
expansions. This model consists of a honeycomb lattice of clusters of four spin-1/2 moments, and 
contains short-range interactions up to six-spin(or eight-spin) terms. The spin in the Kitaev model 
is represented not as these spin-1/2 moments, but as pseudo-spin of the two-dimensional spin singlet 
sector of the four antiferromagnetically coupled spin-1/2 moments within each cluster. Spin corre- 
lations in the Kitaev model are mapped to dimer correlations or spin-chirality correlations in this 
model. This exact construction is quite general and can be used to make other interesting spin-1/2 
models from spin rotation invariant Hamiltonians. We discuss two possible routes to generate the 
high order spin interactions from more natural couplings, which involves perturbative expansions 
thus breaks the exact mapping, although in a controlled manner. 
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I. INTRODUCTION. 



Kitaev's exactly solvable spin-1/2 honeycomb lattice 
model^ (noted as the Kitaev model hereafter) has in- 
spired great interest since its debut, due to its exact 
solvability, fractionalized excitations, and the potential 



to realize non-Abelian anyons. The model simply reads 

a; — links <jk> y — links <jk> 

2— links <jk> 

(1) 

where r^^y^^ are Pauli matrices, and x, y, z-links are de- 
fined in FIG. [H It was shown by Kitaev^ that this spin- 
1/2 model can be mapped to a model with one Majo- 
rana fermion per site coupled to Ising gauge fields on the 
links. And as the Ising gauge fiux has no fluctuation, the 
model can be regarded as, under each gauge fiux config- 
uration, a free Majorana fermion problem. The ground 
state is achieved in the sector of zero gauge fiux through 
each hexagon. The Majorana fermions in this sector have 
Dirac-like gapless dispersion resembling that of graphene, 
as long as \ Jx\i and | J^j satisfy the triangular rela- 
tion, sum of any two of them is greater than the third 
one^. It was further proposed by Kitaev'^ that opening of 
fermion gap by magnetic field can give the Ising vortices 
non-Abelian anyonic statistics, because the Ising vortex 
will carry a zero-energy Majorana mode, although mag- 
netic field destroys the exact solvability. 

Great efforts have been invested to better understand 
the properties of the Kitaev model. For example, sev- 
eral groups have pointed out that the fractionalized Ma- 
jorana fermion excitations may be understood from the 
more familiar Jordan- Wigner transformation of ID spin 
systems^ii^. The analogy between the non-Abelian Ising 
vortices and vortices in p + ip superconductors has been 
raised in serveral works^^— . Exact diagonalization has 
been used to study the Kitaev model on small lattices^. 
And perturbative expansion methods have been devel- 
oped to study the gapped phases of the Kitaev-type 
models^. 

Many generalizations of the Kitaev model have been 
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FIG. 1: The honeycomb lattice for the Kitaev model. Filled 
and open circles indicate two sublattices. x, y, z label the links 
along three different directions used in U]). 

derived as well. There have been several proposals to 
open the fermion gap for the non-Abelian phase without 
spoiling exact solvability^*^. And many generahzations 
to other(even 3D) lattices have been developed in the 
last few year&i2,TiS. AH these efforts have significantly 
enriched our knowledge of exactly solvable models and 
quantum phases of matter. 

However, in the original Kitaev model and its later 
generalizations in the form of spin models, spin rotation 
symmetry is explicitly broken. This makes them harder 
to realize in solid state systems. There are many pro- 
posals to realized the Kitaev model in more controllable 
situations, e.g. in cold atom optical latticesilii^, or in 
superconducting circuit&i^. But it is still desirable for 
theoretical curiosity and practical purposes to realize the 
Kitaev-type models in spin rotation invariant systems. 

In this paper we realize the Kitaev honeycomb lattice 
model as the low energy Hamiltonian for a spin rotation 
invariant system. The trick is not to use the physical spin 
as the spin in the Kitaev model, instead the spin-1/2 in 
Kitaev model is from some emergent two-fold degener- 
ate low energy states in the elementary unit of physical 
system. This type of idea has been explored recently by 
Jackeli and Khaliullin— , in which the spin-1/2 in the Ki- 
taev model is the low energy Kramers doublet created by 
strong spin-orbit coupling of t2g orbitals. In the model 
presented below, the Hilbert space of spin-1/2 in the Ki- 
taev model is actually the two dimensional spin singlet 
sector of four antiferromagnetically coupled spin-1/2 mo- 
ments, and the role of spin-1/2 operators(Pauli matrices) 
in the Kitaev model is replaced by certain combinations 
of Sj • Sfe [or the spin-chirality Sj • (S/j x )] between the 
four spins. 

One major drawback of the model to be presented is 
that it contains high order spin interactions(involves up 
to six or eight spins), thus is still unnatural. However it 
opens the possibility to realize exotic (exactly solvable) 
models from spin-1/2 Hamiltonian with spin rotation in- 
variant interactions. We will discuss two possible routes 
to reduce this artificialness through controlled perturba- 
tive expansions, by coupling to optical phonons or by 
magnetic couplings between the elementary units. 

The outline of this paper is as follows. In Section [H] 
we will lay out the pseudo-spin-1/2 construction. In Sec- 




FIG. 2: Left: the physical spin lattice for the model ((8]). The 
dash circles are honeycomb lattice sites, each of which is ac- 
tually a cluster of four physical spins. The dash straight lines 
are honeycomb lattice bonds, with their type x, y, z labeled. 
The interaction between clusters connected by x, y, z bonds 
are the Jx,y,z terms in ([8]) or ((9]) respectively. Note this is not 
the 3-12 lattice used in Ref.^'^''. Right: enlarged picture of 
the clusters with the four physical spins labeled as 1, ... ,4. 
Thick solid bonds within one cluster have large antiferromag- 
netic Heisenberg coupling Jciustor- 

tion mil the Kitaev model will be explicitly constructed 
using this formalism, and some properties of this con- 
struction will be discussed. In Section|IV]we will discuss 
two possible ways to generate the high order spin in- 
teractions involved in the construction of Section IIIII by 
perturbative expansions. Conclusions and outlook will 
be summarized in Section IVl 

II. FORMULATION OF THE PSEUDO-SPIN-1/2 
FROM FOUR-SPIN CLUSTER. 

In this Section we will construct the pseudo-spin-1/2 
from a cluster of four physical spins, and map the phys- 
ical spin operators to pseudo-spin operators. The map- 
ping constructed here will be used in later Sections to 
construct the effective Kitaev model. In this Section we 
will work entirely within the four-spin cluster, all unspec- 
ified physical spin subscripts take values 1, . . . , 4. 

Consider a cluster of four spin-1/2 moments(called 
physical spins hereafter), labeled by Si_..._4, antiferro- 
magnetically coupled to each other (see the right bot- 
tom part of FIG. [J). The Hamiltonian within the clus- 
ter(up to a constant) is simply the Heisenberg antiferro- 
magnetic(AFM) interactions, 

^^clustor = (^clustcr/2) (Si + S2 + S3 + S4)^ (2) 

The energy levels should be apparent from this form: 
one group of spin-2 quintets with energy SJciustor, three 
groups of spin-1 triplets with energy Jciustcr, and two spin 
singlets with energy zero. We will consider large positive 
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"/cluster limit. So only the singlet sector remains in low 
energy. 

The singlet sector is then treated as a pseudo-spin- 1/2 
Hilbert space. From now on we denote the pseudo-spin- 
1/2 operators as T = (l/2)r, with r the Pauli matri- 
ces. It is convenient to choose the following basis of the 
pseudo-spin 

V- = ±1) = ^(i im)+^-^'\ ;nt) +^"1 ;m> 
+ 1 m;> + c^^'l nn> + mt) 

(3) 



physical spin 


pseudo-spin 


Pi2 , and P34 




Pi3, and P24 


-(l/2)r^ + (V3/2)r^ 


Pi4, and P23 


-(l/2)r^ - (^/3/2)r^ 


-X234, X341, -X412, and X123 


(V3/4)r^ 



TABLE L Correspondence between physical spin operators 
and pseudo-spin operators in the physical spin singlet sector of 
the four antiferromagnetically coupled physical spins. Pjk ~ 
2Sj • Sfc -I- 1/2 are permutation operators, Xjkt ~ Sj • (Sfe x Se) 
are spin-chirality operators. Note that several physical spin 
operators may correspond to the same pseudo-spin operator. 



where w = e^'^'/^ is the complex cubic root of unity, 
I iitt) a-iid other states on the right-hand-side(RHS) are 
basis states of the four-spin system, in terms of quan- 
tum numbers of physical spins 1 , . . . , 4 in sequential or- 
der. This pseudo-spin representation has been used by 
Harris et al. to study magnetic ordering in pyrochlore 
antiferromagnets^ ^ . 

We now consider the effect of Heisenberg-type inter- 
actions Sj ■ Sfc inside the physical singlet sector. Note 
that since any Sj • Sfc within the cluster commutes with 
the cluster Hamiltonian TJciuster @, their action do not 
mix physical spin singlet states with states of other total 
physical spin. This property is also true for the spin- 
chirality operator used later. So the pseudo-spin Hamil- 
tonian constructed below will be exact low energy Hamil- 
tonian, without truncation errors in typical perturbation 
series expansions. 

It is simpler to consider the permutation operators 
Pjk = 2Sj ■ Sfc + 1/2, which just exchange the states 
of the two physical spin- 1/2 moments j and k (j ^ k). 
As an example we consider the action of P34, 



However there is another simpler representation of r^, 
by the spin-chirality operator Xjki = Sj • (S^ x S^). Ex- 
plicit calculation shows that the effect of S2 • (S3 x S4) is 
— (\/3/4)t^ in the physical singlet sector. This can also 
be proved by using the commutation relation [S2 • S3, S2 • 
S4] — iS2 ■ (S3 x S4). A complete list of all chirality 
operators is given in TABLE H] Therefore we can choose 
another representation of , 

= -X234/(V3/4) = -(4/V3)S2 • (S3 x S4) (6) 

The above representations of r^'^'^ are all invariant under 
global spin rotation of the physical spins. 

With the machinery of equations (O, and ([5]), it 
will be straightforward to construct various pseudo-spin- 
1/2 Hamiltonians on various lattices, of the Kitaev vari- 
ety and beyond, as the exact low energy effective Hamil- 
tonian of certain spin- 1/2 models with spin- rotation sym- 
metry. In these constructions a pseudo-spin lattice site 
actually represents a cluster of four spin-1/2 moments. 



P3A\r' = -1) = ^(1 ;;n> ;m> +^^'1 ;nt) 
+ itm>+'^imt>+'^'inn> 

= = +1) 



and similarly P34|t^ = —1) = = -fl). Therefore P34 
is just in the physical singlet sector. A complete list 
of all permutation operators is given in TABLE HI We 
can choose the following representation of and r^. 



= P12 - 2Si • S2 + 1/2 
ry = (Fi3 - Pi4)/V3 - (2/%/3)Si • (S3 - S4) 



(4) 



Many other representations are possible as well, because 
several physical spin interactions may correspond to the 
same pseudo-spin interaction in the physical singlet sec- 
tor, and we will take advantage of this later. 

For we can use = — ir^r*', where i is the imagi- 
nary unit, 

= -i(2/V3)(2Si • S2 + l/2)Si • (S3 - S4) (5) 



III. REALIZATION OF THE KITAEV MODEL. 

In this Section we will use directly the results of the 
previous Section to write down a Hamiltonian whose low 
energy sector is described by the Kitaev model. The 
Hamiltonian will be constructed on the physical spin lat- 
tice illustrated in FIG. [2] In this Section we will use 
J, k to label four-spin clusters (pseudo-spin- 1/2 sites), the 
physical spins in cluster j are labeled as S^i, . . . , Sj4. 

Apply the mappings developed in Section [Hi we have 
the desired Hamiltonian in short notation, 



-^cluster 



cluster 



x — links <jk> 

E 

y— links <jk> 2; — links <jk> 



E 



1 T^T^ 



(7) 



where j, k label the honeycomb lattice sites thus the four- 
spin clusters, i?ciustcr is given by ([2]), r^'?''^ should be 
replaced by the corresponding physical spin operators in 
(HI) and ([5]) or ([H]), or some other equivalent representa- 
tions of personal preference. 
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Plug in the expressions ^ and ([6]) into (O , the Hamil- tonian reads explicitly as 
I 



H = 5](Jclustcr/2)(Sji + S,2 + S,3 + S,4)' ' (16/9) [S,2 • (8,3 X S,4)][Sfe2 ' (Sfc3 X 8^4)] 

j 2 — links <jk~> 

- ^(2S,i-S,2 + l/2)(2Sfei.Sfe2 + l/2)- ^ J,(4/3)[S,i.(S,3-S,4)][Sfei-(Sfc3-Sfc4)] 

X— links <jk> y — links <jk> 

(8) 



r 



While by the represenation (j4]) and the Hamilto- nian becomes 
I 



luster /2)(S,i+S,2 + S,3+S,4)' 

j 

J,(2S,i-S,2 + l/2)(2Sfci-Sfe2 + l/2)- (4/3) [S,i .(8,3-8^4)] [Sfci-(Sfc3-Sfc4)] 

x — links <jk> y — links <jk> 

- ^ J, (-4/3)(2S,3 • S,4 + l/2)[S,i • (S,3 - S,4)](2Sfe3 • 8^4 + l/2)[Sfei • (Sfc3 - 8^4)] 

2 — links <jk> 

(9) 



r 



This model, in terms of physical spins S, has full 
spin rotation symmetry and time-reversal symmetry. A 
pseudo-magnetic field term h ■ Tj term can also be 
included under this mapping, however the resulting Ki- 
taev model with magnetic field is not exactly solvable. 
It is quite curious that such a formidably looking Hamil- 
tonian ([8|), with biquadratic and six-spin(or eight-spin) 
terms, has an exactly solvable low energy sector. 

We emphasize that because the first intra-clustcr term 
X^ciustcr -^cluster commutcs with the latter Kitaev terms 
independent of the representation used, the Kitaev model 
is realized as the exact low energy Hamiltonian of this 
model without truncation errors of perturbation theories, 
namely no (| Jx,j,,z|/ ^ciustor)^ or higher order terms will 
be generated under the projection to low energy clus- 
ter singlet space. This is unlike, for example, the t/U 
expansion of the half-filled Hubbard mode P^i^^ , where 
at lowest t'^/U order the effective Hamiltonian is the 
Heisenberg model, but higher order terms {t^ /U^ etc.) 
should in principle still be included in the low energy ef- 
fective Hamiltonian for any finite t/U . Similar compari- 
son can be made to the perturbative expansion studies of 
the Kitaev-type models by Vidal et al.^, where the low 
energy effective Hamiltonians were obtained in certian 
anisotropic (strong bond/triangle) limits. Although the 
spirit of this work, namely projection to low energy sec- 
tor, is the same as all previous perturbative approaches 
to effective Hamiltonians. 



Note that the original Kitaev model ([T]) has three- 
fold rotation symmetry around a honeycomb lattice site, 
combined with a three-fold rotation in pseudo-spin space 
(cyclic permutation of r^, r^, r^). This is not apparent 
in our model ([S]) in terms of physical spins, under the 
current representation of r^'^'^. We can remedy this by 
using a different set of pseudo-spin Pauli matrices t'^'^'^ 
in dll), 




With proper representation choice, they have a symmet- 
ric form in terms of physical spins. 



r'y 




-(4/3)82 

-(4/3)83- 

-(4/3)84- 



So the symmetry mentioned above can be realized by a 
three- fold rotation of the honeycomb lattice, with a cyclic 
permutation of 82, 83 and 84 in each cluster. This is in 
fact the three-fold rotation symmetry of the physical spin 
lattice illustrated in FIG. [2] However this more symmet- 
ric representation will not be used in later part of this 
paper. 
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Another note to take is that it is not necessary to have 
such a highly symmetric cluster Hamiltonian ([2]). The 
mappings to pseudo-spin- 1/2 should work as long as the 
ground states of the cluster Hamiltonian are the two-fold 
degenerate singlets. One generalization, which conforms 
the symmetry of the lattice in FIG. [21 is to have 



cluster 



(>/clustor/2)(r • Si + S2 + S3 



(11) 



with Jciustcr > and < r < 3. However this is not 
convenient for later discussions and will not be used. 

We briefly describe some of the properties of (|8]). Its 
low energy states are entirely in the space that each of the 
clusters is a physical spin singlet (called cluster singlet 
subspace hereafter) . Therefore physical spin correlations 
are strictly confined within each cluster. The excitations 
carrying physical spin are gapped, and their dynamics 
are 'trivial' in the sense that they do not move from one 
cluster to another. But there are non-trivial low energy 
physical spin singlet excitations, described by the pseudo- 
spins defined above. The correlations of the pseudo-spins 
can be mapped to correlations of their corresponding 
physical spin observables (the inverse mappings are not 
unique, c.f. TABLE |l|. For example r^'^ correlations 
become certain dimer-dimer correlations, correlation 
becomes chirality-chirality correlation, or four-dimer cor- 
relation. It will be interesting to see the corresponding 
picture of the exotic excitations in the Kitaev model, e.g. 
the Majorana fermion and the Ising vortex. However this 
will be deferred to future studies. 

It is tempting to call this as an exactly solved spin liq- 
uid with spin gap (^ Jciustcr), an extremely short-range 
resonating valence bond(RVB) state, from a model with 
spin rotation and time reversal symmetry. However it 
should be noted that the unit cell of this model contains 
an even number of spin-1/2 moments (so does the orig- 
inal Kitaev model) which does not satisfy the stringent 
definition of spin liquid requiring odd number of elec- 
trons per unit cell. Several parent Hamiltonians of spin 
liquids have already been constructed. See for example, 
Ref^-^. 



IV. GENERATE THE HIGH ORDER PHYSICAL 
SPIN INTERACTIONS BY PERTURBATIVE 
EXPANSION. 

One major drawback of the present construction is that 
it involves high order interactions of physical spins [see 
(|S]) and (|n])], thus is 'unnatural'. In this Section we will 
make compromises between exact solvability and natu- 
ralness. We consider two clusters j and k and try to 
generate the Jx,y,z interactions in ([T]) from perturbation 
series expansion of more natural(lower order) physical 
spin interactions. Two different approaches for this pur- 
pose will be laid out in the following two Subsections. In 
Subsection IIV Al we will consider the two clusters as two 
tetrahedra, and couple the spin system to certain opti- 
cal phonons, further coupling between the phonon modes 
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FIG. 3: Illustration of the tetragonal to orthorhombic 
Qi'(top) and (bottom) distortion modes, (a) Perspective 
view of the tetrahedron. 1, . . . , 4 label the spins. Arrows in- 
dicate the motion of each spin under the distortion mode, (b) 
Top view of (a). (c)(d) Side view of (a). 



of the two clusters can generate at lowest order the de- 
sired high order spin interactions. In Subsection IIV Bl we 
will introduce certain magnetic, e.g. Heisenberg-type, in- 
teractions between physical spins of different clusters, at 
lowest order (second order) of perturbation theory the de- 
sired high order spin interactions can be achieved. These 
approaches involve truncation errors in the perturbation 
series, thus the mapping to low energy effect Hamilto- 
nian will no longer be exact. However the error intro- 
duced may be controlled by small expansion parameters. 
In this Section we denote the physical spins on cluster 
jik) as jT, . . . , j4 (fcl, . . . , fc4), and denote pseudo-spins 
on cluster j(k) as Tj(T}.). 



A. Generate the High Order Terms by Coupling to 
Optical Phonon. 

In this Subsection we regard each four-spin cluster 
as a tetrahedron, and consider possible optical phonon 
modes(distortions) and their couplings to the spin sys- 
tem. The basic idea is that the intra-cluster Heisen- 
berg coupling Jciuster can linearly depend on the dis- 
tance between physical spins. Therefore certain distor- 
tions of the tetrahedron couple to certain linear combi- 
nations of S^ • Sm- Integrating out phonon modes will 
then generate high order spin interactions. This idea has 
been extensively studied and applied to several magnetic 
materials^""—. More details can be found in a recent 
review by Tchernyshyov and Chern^^. And we will fre- 
quently use their notations. In this Subsection we will 
use the representation ([5]) for t^. 

Consider first a single tetrahedron with four spins 
1, . . . , 4. The general distortions of this tetrahedron can 
be classified by their symmetry (see for example Refi^) . 
Only two tetragonal to orthorhombic distortion modes, 
Qf and Q2 (illustrated in FIG. [3]), couple to the pseudo- 
spins defined in Section [Hi A complete analysis of all 
modes is given in Appendix \^ The coupling is of the 
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form 



where J' is the derivative of Heisenberg couphng Jciustcr 
between two spins i and m with respect to their distance 
rim, J' = dJciustor/dr^m; are the generahzed coor- 
dinates of these two modes; and the functions ff'2 are 



fi 



= (1/2)(S2 • S4 + Si • S3 - Si • S4 - S2 • S3), 
- 71712(81 • S4 + S2 • S3 + S2 • S4 + Si • S3 
— 2Si • S2 ^ 2S3 • S4). 



According to TABLE U we have /f = -{V^/^^ and 
f2 — ("\/3/2)t^. Then the couphng becomes 



(V3/2)J'(-Qfr-+Qfr^) 



(12) 



The spin-lattice(SL) Hamihonian on a single cluster j 
is [equation (1.8) in Ref.— ], 



-^cluster J, SL — -^cluster ^ ^~ 



k 
2 



(Qf,)^ 



(13) 



where A: > is the elastic constant for these phonon 
modes, J' is the spin-lattice coupling constant, Qfj and 
Q2j are the generalized coordinates of the Qf and 
distortion modes of cluster j, ffciustcr j is ([2]). As al- 
ready noted in Refi^, this model does not really break 
the pseudo-spin rotation symmetry of a single cluster. 

Now we put two clusters j and k together, and in- 
clude a perturbation A ^perturbation to the optical phonon 
Hamiltonian, 

Hjk,SL =^^clustcr j, SL + -f^clustor fc, SL 

+ A -ffpcrturbation[Qij, Q2j' Qlfe' Q2k\ 



where A (in fact A/fc) is the expansion parameter. 

Consider the perturbation TJperturbation = Qfj ■ Qf^, 
which means a coupling between the Qf distortion 
modes of the two tetrahedra. Integrate out the optical 
phonons, at lowest non-trivial order, it produces a term 
{3 J''^ X)/{Ak'^)T^ ■ rf. This can be seen by minimizing 
separately the two cluster Hamiltonians with respect to 
Qf, which gives Qf = (^3 J')/(2 fc)r^, then plug this 
into the perturbation term. Thus we have produced the 
Jx term in the Kitaev model with ~ —(3 J'^ A)/(4 k^). 

Similarly the perturbation i/pcrturbation = Qfj ' Qik 
will generate (3 J'^ A)/(4 k'^) tJ ■ at lowest non-trivial 
order. So we can make Jy — — (3 J''^ A)/(4 fc^). 

The ■ t| coupling is more difficult to get. We 

By the above reasoning, we 



J J 



•k'k- 



3 

treat it as 

J J n. (V 

need an anharmonic coupling i/pcrturbation = QijQfj 
QfkQik- It will produce at lowest non-trivial or- 



der (9 J'4A)/(16fc4)r;rJ 



'k 'k- 



Thus we have 



(9 J'^\)/{16k^). 

Finally we have made up a spin-lattice model ifsL, 
which involves only S^ • S^ interaction for physical spins. 



cluster 



HsL — ^ -^cluster, SL + ^ Qlj ' Q 

2; — links <jk> 

Ay Q2J ■ Q^ 

y — links <jk> 

QijQ2j ■ QikQ 

2— links <jk> 



l2k 



i2k 



where Qf^ is the generalized coordinate for the Qf mode 



on cluster j, and Qff^, Qij' Qik ^^'^ similarly defined; 



y = -{4:Jx^yk^)/{3J'^) and A^ = {16J,k^)/i9J'^); the 
single cluster spin-lattice Hamiltonian i/ciustcr, sl is (|13p . 

Collect the results above we have the spin-lattice 
Hamiltonian iJgL explicitly written as, 



i^SL- [(•^clustcr/2)(S,i+S,2 + S,3+Sj4)2 + ^(Qf^.)' + f(0fj) 

'■( 



cluster j 



^ V12 

nE ■ Sj4 + Sji • Sj3 — Sji • Sj4 — Sj2 • Sj3 ^ ^ ^ ^ 



2 



3J/2 ■ Qlk - 2^ 3J,2 ■ Q2k + 2^ gj/4 *5ljQ2j ' QlkQ2k 

X— links <jk> y — links <jk> 2 — links <jk> 



The single cluster spin-lattice Hamiltonian [first three harmonic (on x- and y-links of honeycomb lattice) and an- 
lines in (|14p ] is quite natural. However we need some harmonic coupling (on z-links) between optical phonon 
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modes of neighboring tetrahedra. And these couphng 
constants Xx,y,z need to be tuned to produce Jx,y,z of 
the Kitaev modeL This is stih not easy to implement in 
sohd state systems. At lowest non-trivial order of pertur- 
bative expansion, we do get our model Higher order 
terms in expansion destroy the exact solvability, but may 
be controlled by the small parameters \x,y,zlk. 

B. Generate the High Order Terms by Magnetic 
Interactions between Clusters. 

In this Subsection we consider more conventional per- 
turbations, magnetic interactions between the clusters, 
e.g. the Heisenberg coupling Sj • with j and k belong 
to different tetrahedra. This has the advantage over the 
previous phonon approach for not introducing additional 
degrees of freedom. But it also has a significant disad- 
vantage: the perturbation does not commute with the 
cluster Heisenberg Hamiltonian ©, so the cluster sin- 
glet subspace will be mixed with other total spin states. 
In this Subsection we will use the spin-chirality represen- 
tation ^ for . 

Again consider two clusters j and k. For simplicity 
of notations define a projection operator Vjk = 'Pj'Pk, 
where Vj^k is projection into the singlet subspace of clus- 
ter j and k, respectively, Vj^k = Es=±i \'^lk = ■^X'^/.fe = 
s\. For a given perturbation A Ifporturbation with small 
parameter A (in factor A/ Jdustcr is the expansion param- 
eter), lowest two orders of the perturbation series are 

A T^j/c -/^pcrturbation'T^j/i.' ^" A 7^j/c-/^pcrturbation(l T^jk) 
X [0 -^^clustcr j -^cluster k] (1 T^jk^-^pciturha,tionT^jk 

(15) 

With proper choice of A and iJpcrturbation we can generate 



the desired Jx,y,z terms in ([8]) from the first and second 
order of perturbations. 

The calculation can be dramatically simplified by the 
following fact that any physical spin-1/2 operator 5^'^'^ 
converts the cluster spin singlet states |t^ = ±1) into 
spin-1 states of the cluster. This can be checked by 
explicit calculations and will not be proved here. For 
all the perturbations to be considered later, the above 
mentioned fact can be exploited to replace the factor 
[0 - i?ciustor J ~ ^^ciustcr k]~^ in the second order pertur- 
bation to a c-number (—2 Jdustcr) 

The detailed calculations are given in Appendix |B] We 
will only list the results here. 

The perturbation on x-links is given by 

Aa; -^perturbation, x — Xx[Sji ■ Sfci +sgn{Jx) ■ (Sj2 • Sfe2)] 

— Jx{Sji ■ Sj2 + Ski ■ Sfe2)- 



where A^; = ^12\Jx\- Jdustcr, sgn(Jj;) = ±1 is the sign 
of Jx- 

The perturbation on y-links is 

Ay J^pcrturbation, y 

^Xy[Sji ■ Sfci + sgn(Jj^) • (Sj3 - Sji) ■ (Sfc3 - Sm)] 

~ I'A/I(^j3 • Sj4 -I- Sfe3 • Sfc4) 



with Xy = \/4:\Jy \ ■ Jdustcr. 

The perturbation on z-links is 

Xz -^^perturbation, z 
= Az[Sj2 • (Sfe3 X Sfe4) + Sgn{J^) ■ Sk2 ■ (Sj3 X Sji) 

— \Jz\{Sj3 ■ Sj4 -I- Sfe3 • Sm)- 



with Xz = 4a/| J^l • Jdustcr- 

The entire Hamiltonian iJmagnctic reads explicitly as, 



-^magnetic — / ^ (>/clustcr/2) (S ji + S j2 + Sj3 + Sj4 )^ 
cluster j 

+ ^ {y/l2\Jx\ ■ Jduster[Sjl • Sfci + Sgn(J^) - {Sj2 ■ 8^2)] - Jx{Sjl ■ Sj2 + Sfci • Sfe2)} 

a; — links <jk> 

+ ^ {y^4|Jy| ■ Jdustcr [Sji • (S/e3 — S^) + Sgn(Jy)S/el * (Sj3 — 8^4)] — |Jy|(Sj3 • Sj4 + 8/^3 ' ^kA.)} 

y— links <jk> 

+ ^ {4v/|J2| ■ Jdustcr [Sj2 • (Sfc3 X Sfc4) + Sgn(J2)Sfc2 • (Sjs X S^^)] - \ Jz\iS.j3 ■ Sj4 + Sfc~3 • Sfe4)}. 

links <jk> 

I 



In we have been able to reduce the four spin in- 
teractions in ^ to inter-cluster Heisenberg interactions, 
and the six-spin interactions in to inter-cluster spin- 
chirality interactions. The inter-cluster Heisenberg cou- 
plings in i/pcrturbation x,y may be easier to arrange. The 



inter-cluster spin-chirality coupling in iJpcrturbation z ex- 
plicitly breaks time reversal symmetry and is probably 
harder to implement in solid state systems. However 
spin-chirality order may have important consequences 
in frustrated magnets^i^, and a realization of spin- 



8 



chirality interactions in cold atom optical lattices has 
been proposectf^. 

Our model ([S]) is achieved at second order of the per- 
turbation series. Higher order terms become trunca- 
tion errors but may be controlled by small parameters 

^x^y^z / '-^cluster ^ \l \ Jx,y ,z I / "-^cluster ■ 

V. CONCLUSIONS. 

We constructed the exactly solvable Kitaev honeycomb 
modeli as the exact low energy effective Hamiltonian of 
a spin-1/2 model [equations ([5]) or with spin-rotation 
and time reversal symmetry. The spin in Kitaev model is 
represented as the pseudo-spin in the two-fold degenerate 
spin singlet subspace of a cluster of four antiferromag- 
netically coupled spin-1/2 moments. The physical spin 
model is a honeycomb lattice of such four-spin clusters, 
with certain inter-cluster interactions. The machinery 
for the exact mapping to pseudo-spin Hamiltonian was 
developed (see e.g. TABLE which is quite general 
and can be used to construct other interesting (exactly 
solvable) spin-1/2 models from spin rotation invariant 
systems. 

In this construction the pseudo-spin correlations in the 
Kitaev model will be mapped to dimer or spin-chirality 
correlations in the physical spin system. The correspond- 
ing picture of the fractionalized Majorana fermion exci- 
tations and Ising vortices still remain to be clarified. 

This exact construction contains high order physical 
spin interactions, which is undesirable for practical im- 
plementation. We described two possible approaches to 
reduce this problem: generating the high order spin in- 
teractions by perturbative expansion of the coupling to 
optical phonon, or the magnetic coupling between clus- 
ters. This perturbative construction will introduce trun- 
cation error of perturbation series, which may be con- 
trolled by small expansion parameters. Whether these 
constructions can be experimentally engineered is how- 
ever beyond the scope of this study. It is conceivable that 
other perturbative expansion can also generate these high 
order spin interactions, but this possibility will be left for 
future works. 
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Appendix A: Coupling between Distortions of a 
Tetrahedron and the Pseudo-spins 

In this Appendix we reproduce from Ref.^^ the cou- 
plings of all tetrahedron distortion modes to the spin 



system. And convert them to pseudo-spin notation in 
the physical spin singlet sector. 

Consider a general small distortion of the tetrahedron, 
the spin Hamiltonian becomes 

-^cluster, SL — ("^cluster 

i i<m 

(Al) 

where Srgm is the change of bond length between spins 
£ and to, J' is the derivative of Jdustcr with respect to 
bond length. 

There are six orthogonal distortion modes of the tetra- 
hedron [TABLE 1.1 in Ref.35]. One of the modes A is the 
trivial representation of the tetrahedral group Td] two E 
modes form the two dimensional irreducible representa- 
tion of Td] and three T2 modes form the three dimen- 
sional irreducible representation. The E modes are also 
illustrated in FIG. |1 

The generic couplings in (|Aip [second term] can be 
converted to couplings to these orthogonal modes, 

J'iQ^f^ + Qf /f + Qifi + Q?f? + Q?f? + Ql'f?) 

where Q are generalized coordinates of the corresponding 
modes, functions / can be read o ff fro m TABLE 1.2 of 
Ref.35. For the A mode, Srg^ = y^Q^, so is 

= 7273(81 • S2 + S3 • S4 + Si • S3 

-|- S2 ■ S4 + Si • S4 + S2 ■ S3). 

The functions f^2 the E modes have been given before 
but are reproduced here, 

fi = (1/2)(S2 • S4 + Si ■ S3 - Si • S4 - S2 • S3), 
/f = v/l7l2(Si • S4 + S2 • S3 + S2 • S4 + Si • S3 
— 2Si • S2 — 2S3 • S4). 

The functions 3 fo^' the T2 modes are 





= (S2 


•S3 


-Si 


•S4 


J2 


-(Si 


•S3 


-S2 


•S4 


J3 


= (Si 


•S2 


-S3 


•S4 



Now we can use TABLE |I] to convert the above cou- 
plings into pseudo-spin. It is easy to see that and 
/?2 3 when converted to pseudo-spins, namely 

projected to the physical spin singlet sector. But fi — 

(Pl4 + P23 + P24 + Pl3-2Pl2-2P34)/(4V3) = -{V^/2)t^ 

and fi = (P24 + Pi3 - Pi4 - P23)/4 = {V5/2)Ty. This 
has already been noted by Tchernyshyov et al.^-, only 
the E modes can lift the degeneracy of the physical spin 
singlet ground states of the tetrahedron. Therefore the 
general spin lattice coupling is the form of (jl2p given in 
the main text. 
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Appendix B: Derivation of the Terms Generated by 
Second Order Perturbation of Inter-cluster 
Magnetic Interactions 

In this Appendix we derive the second order pertur- 
bations of inter-cluster Heisenberg and spin-chiraUty in- 
teractions. The results can then be used to construct 

(USD. 

First consider the perturbation A ^perturbation = H^ji ' 
Sfci + ?'(Sj2 ■ Sfe2)]j where r is a real number to be tuned 
later. Due to the fact mentioned in Subsection IIVB| 
the action of iJperturbation On any cluster singlet state 
will produce a state with total spin-1 for both cluster j 
and k. Thus the first order perturbation in (fTSj) van- 
ishes. And the second order perturbation term can be 
greatly simplified: operator (1 — 'Pjk)[0 — ^^ciustor j — 

-^^clustcr k 

] ^(1 - Vjk) can be replaced by a c-number 
(— 2Jciustor)~^- Therefore the perturbation up to second 
order is 



A2 



2J< 



cluster 



" (-^perturbation) ^jk 



This is true for other perturbations considered later in 
this Appendix. The cluster j and cluster k parts can be 
separated, this term then becomes (a, b = x^y, z), 



A2 



2J< 



cluster 



E 

a.b 



[P,S'^,S],VyVkSi,Sl,Vk 



+ 2rV,S'^,S]^VrVkSl,Sl^Vk 
+ r^VjS';^S]^VyVkSt^Sl^Vk] 

Then use the fact that VjSf^S^^Vj = Sab{l/i)'Pj{Sjt ■ 
Sjm)'Pj by spin rotation symmetry, the perturbation be- 
comes 



A2 



6</clustcr 

A2 



9 + 9r2 



6Jc 



luster 



16 
9 + 9r^ 



16 



+ 2rVjkiSji ■ Sj2)(Sfei • Sk2)Vjk 
+ {r/2)TfT^ r/2 



rVjki^ji ■ Sj2 + Ski ■ Sk2)Vjk 



So we can choose — (r A^)/(12 Jduster) — ~Jx, and include 
the last intra-cluster Sji • Sj2 + S^i • 8^2 term in the first 
order perturbation. 

The perturbation on x- links is then (not unique), 

Aa; -^^perturbation, x 

=Ax[Sji • Sfci + sgn(Ja;) • (Sj2 • Sfc2)] 

— Jx{Sjl ■ Sj2 + Ski ■ Sfe2) 



with Xj; = a/12| Ja;| • Jciuster, and r = sgn(Jj;) is the sign 
of Jx- The non-trivial terms produced by up to second 
order perturbation will be the TjT^ term. Note that the 
last term in the above equation commutes with cluster 
Hamiltonians so it does not produce second or higher 
order perturbations. 



Similarly considering the following perturbation on y- 

links, A -ffpcrturbation = ^[Sji • (Sfc3 — Ski) +rSkl ' (Sj3 — 

Sj4)]. Following similar procedures we get the second 
order perturbation from this term 



A2 



9r^ 



6t/clustc 

+ 2rVjk[Sji-{S^3- 
— {3/2)VjkiSk3 ■ Sk4 
A2 r9 + 9r2 



Sj4)][SA:i • (Sfc3 — Ski)]^^ 



jk 



Sj3 • Sj/CjVjk 



6J, 



cluster 
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2r (3/4)' 



(3/2) V,k{Sk3 ■ Ski + Sj3 • S,i)V,k 



So we can choose — (r A^)/(4Jciuster) = ^Jyt and include 
the last intra-cluster 8^:3 • SkA + Sj^ ■ Sj^ term in the 
first order perturbation. 

Therefore we can choose the following perturbation on 
y-links (not unique). 

Ay -^^perturbation, y 

^Xy[Sji ■ Sfel + sgn( Jy) • (Sj3 - Sj4) • (Sfc3 - Sfc4)] 

^ l'^l/l(Sj3 • Sji + Sfe3 • Sfc4) 



with Ay = A/4|Jy| • Jciuster, r = sgn(Jy) is the sign of Jy. 

The tJt| term is again more difficult to get. We use 
the representation of by spin-chirality ©. And con- 
sider the following perturbation 

-ffperturbation = Sj2 ■ (Sj3 X 8^4) + r Sfe2 ' (Sj3 X 8^4) 

The first order term in (1151) vanishes due to the same 



reason as before. There are four terms in the second 
order perturbation. The first one is 

\^VjkSj2-{Sk3y-SkA){l-Vjk) 

X [0 -ffclustor j -^cluster fe] 

X (1 - V.jk)Sj2 ■ (Sfe3 X Ski)Vjk 

For the cluster j part we can use the same arguments 
as before, the -ffciuster j can be replaced by a c-number 
"/cluster- For the cluster k part, consider the fact that 
Sfc3 X Sfe4 equals to the commutator — i[Sfc4,Sfe3 • 8^4], 
the action of Sfe3 x 8^4 on physical singlet states of k will 
also only produce spin-1 state. So we can replace the 
/^cluster fe in the denominator by a c-number Jciuster as 
well. Use spin rotation symmetry to separate the j and 
k parts, this term simplifies to 

A^ 

'^j^j^ ' ^j2^j ■ 'Pfe(Sfe3 X Sfc4) • (Sfc3 X SkijVk- 

t> "^cluster 

Use (S)2 = 3/4 and 

{Sk3 X Sfc4) • (S/j3 X Sfe4) 



E 



/ Qa oh oa oh oa ob oh oa \ 
Wk3'-^k4'-*k3'-*k4 '-^k3'^k4'^k3'-^k4) 



{^k3 • Sfc3)(S/,4 ■ S/,4) - ^^S^^Sl^lSahf^ - Sl^S\ 
a,h 

9/16 +(Sfe3.Sfe4)(Sfc3-Sfe4)- (3/8) 
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this term becomes 
A2 



6 Jclust 
^2 



(3/4)[3/16+(tV2- 1/4)2 



= - (A^)/(32Jeiustcr) • (2 - T^). 

Another second order perturbation term r^A^ 'Pjk^k2 ■ 



(S,3 X S,4)(l - Vjk)[Q - 



luste 



H, 



cluster k\ 



-Hi 



r/). 



'Pjfe)Sfe2 • (Sj3 X Sj4)Vjk can be computed in the similar 
way and gives the result — (r^ A^)/(32Jciustcr) • (2 
For one of the cross term 

rX^VjkSj2-{Sk3 xSk4){l-Vjk) 

X [0 — -^cluster j ^ -^^clustor k] 

X (1 - Vjk)Sk2 ■ (Sj3 X Sji)rjk 

We can use the previous argument for both cluster j and 
k, so (l-PAB)[0-i/ciustcr j-ffciustcr k]'^i^-'Pjk) Can be 
replace by c-number (— 2 Jciustoi)^^- This term becomes 



rA2 



2J( 



cluster 



-PjfciSja • (Sfc3 X SM)][Sfe2 • (S,3 X S,3)]Vjk- 



Spin rotation symmetry again helps to separate the terms 
for cluster j and k, and we get — (r A^)/(32 Jciustcr) • '^/'^fc ■ 
The other cross term r \^ Vjk^k2 ■ (Sj3 x Sj4)(l — 

'Pjk)[0 — -ffclustor j — ^^clustor fe]^ (1 — 'Pjk)^j2 ' (Sfc3 X 

^k4)'Pjk gives the same result. 

In summary the second order perturbation from A[Sj2 • 

(Sj3 X Sj4) + rSfc2 ■ (Sj3 X Sj4)] is 



rA2 



A2 



16Jc 



luste 



32 J( 



cluster 



Using this result we can choose the following pertur- 
bation on z-links, 



A^ -^^perturbation, z 

=Xz[Sj2 ■ (Sfc3 X Sm) + sgn(J2) • Sk2 ■ (Sj3 X Sji) 

~ l'^z|(Sj3 • Sj4 + 8^:3 • Sfe4) 



with Az = 4^1 J2I Jciustcr, r = sgn(J2) is the sign of J^. 
The last term on the right-hand-side is to cancel the non- 
trivial terms (r^ -|- rjf )A^/(32 Jciustcr) from the second 
order perturbation of the first term. Up to second order 
perturbation this will produce — J^tJtI interactions. 

Finally we have been able to reduce the high order 
interactions to at most three spin terms, the Hamiltonian 



magnetic 



IS 



^magnetic — ^ ^ -f^cluster j ~l~ ^ ^ Ax-f^perturbation x 

j a; — links <jk> 

^ ^ -^y-^pcrturbation y 

y — links <jk> 

~t~ ^ ^ '^z -^perturbation z 

2 — links <jk> 



where 7?ciuster j are given by {l2|), Aa;^y^2: -ffporturbation x^y^z 

are given above. Plug in relevant equations we get (jl6p 
m Subsection HVB] 
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